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I. INTRODUCTION 

in 

G\ Forb ound conservative Hamiltonian systems, the main 
O^goal of semiclassical mechanics is to obtain estimates for 
'the quantum eigenvalues and eigenfunctions using clas- 
^s^cal trajectories. Since these eigenvalues and eigenfunc- 
, ^ tions are time independent quantities, it is natural to 
^- expect that only classical trajectories that define time in- 
^variant classical manifolds be important semiclassically. 
^ For example, the number of Planck cells in the phase 
space volume enclosed by the energy shell of energy E 
ives a good semiclassical estimate for the number of 
^^eigenstates with energy less than E . 

For chaotic systems, the only time invariant mani- 
f~jblds other than the energy shell are periodic orbits. 
^/^Consequently periodic orbits are at the center of ef- 
^Sforts to obtain viable semiclassical quantization schemes. 
"^Gutzwiller's trace formula [1] uses periodic orbits to ex- 
^rjlain the oscillations in the density of states about the av- 
'^erage density of states. A similar expression was recently 
'developed by Bogomolny [2] to explain the accentuations 
^and attenuations in the average coordinate probability 
fH density about the classical periodic orbits, the so-called 
(Jscars [3]. 

In this paper we present applications of Gutzwiller's 
trace formula and Bogomolny's formula to the following 
generic Hamiltonian system [4]: 
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with n = 0.1. This anharmonic Hamiltonian system has 
been code named Nelson [5] and has been studied exten- 
sively both classically and semiclassically [4-6]. The slow 
progress in the application of periodic orbit theories to 
physical systems such as the Nelson Hamiltonian is due 
to the daunting task of finding all the periodic orbits with 
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periods r less than r max = 27r/}/e, where e is the energy 
resolution desired. We used a method which guarantees 
that all periodic orbits up to a chosen maximum period 
are found. This allowed us to make stringent tests of 
Gutzwiller's trace formula and of Bogomolny's formula. 

The outline of the paper is as follows. In section II 
we review the semiclassical theory, paying particular at- 
tention to symmetry considerations. We present in sec- 
tion III the method we used to find classical periodic 
orbits of a chaotic system in two spatial dimensions and 
then go on to find the desired orbits for the Nelson Hamil- 
tonian. In section III B we use these periodic orbits to 
calculate the semiclassical oscillations in the density of 
states using Gutzwiller's trace formula for different en- 
ergy smoothing parameters e. We then compare our re- 
sults with the exact quantum oscillations in the density of 
states. Finally we present comparisons of the semiclassi- 
cal theory for the oscillating part of the energy smoothed 
probability density with the exact quantum results. 



II. THEORY 

Consider a bound Hamiltonian system in two spatial 
dimensions of the form 



#(q,p) = y + nq)- 



(2) 



(1) Its spectral probability density A(q, E) is defined as [2,7]: 



A(q,E) = J2HE-E k )\(q\^ k 



(3) 



where E k are the eigenvalues of H and are the cor- 
responding eigenstates. The density of states d(E) is 
obtained by integrating the spectral probability density 
over the spatial coordinates q: 



d(E)= j dq A(q,E) 
= Y J 8(E-E k ). 



(4) 
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A semiclassical formula for A(q, E) is obtained by rewrit- 
ing Eq. (3) as: 

1 f°° 

A(q,£) = — Re/ dt(q\e- iHt / n \q)e iEt / n . (5) 

1"" Jo 

An appropriate semiclassical approximation for the prop- 
agator (q\e~ lHt / n \q) is then used and the resulting time 
integral is evaluated by stationary phase. The dominant 
contribution to this integral comes from "zero length" 
trajectories, those whose propagation time is infinitesi- 
mal; this gives rise to the average spectral probability 
density A avg (q_, E). The remaining contributions to the 
time integral come from classical trajectories which start 
at q and come back to q in a nonzero time t; these tra- 
jectories produce the oscillatory part A osc (q, E) of the 
spectral probability density (the so called "scars"). 



A. Average spectral probability density 

The average spectral probability density A avg (q_, E) is 
given by 

A a „, (q, E) = -±-Re f dt <q|e-'^/ B |q)e" 3t / B (6) 

where r denotes a time much smaller than the period 
of the shortest periodic orbit. The Wigner representa- 
tion [8-10] for the evolution operator is used to evaluate 
A a „«,(q, E): 



A a „«,(q, E) 



(2TThf Jo 



Re / dt dp [e 



-mt/n-i jEt/n 



(7) 



where [. . ] w denotes the Weyl symbol of the enclosed 
operator. For short times [8,10] 



-13H-, 



l--A 1 ( (l ,p,(]) + 0(h 4 )} (8) 



where f3 = it/Ti. We shall neglect terms of 0(h 4 ). The 
Weyl symbol for Hamiltonians of the form (2) is simply 
if(q,p). Also [8] 

^i(q, P, /?) = P 2 W - ^ [(W) 2 + (p • V) V] . (9) 

The p x py term in Eq. (9) does not contribute to 
A avg (q, E). An integration by parts of the terms con- 
taining p 2 x or p 2 gives 



xe KE-H) (1Q) 

with 



^"(q,/?) = f/?W-^(W) 2 . 
The momentum and time integrals then give: 
1 



(11) 



A a „«,(q, E) 



2irh z 



dz Ai(z) 



x 



1^V 2 W - j- 

12 UlW|) 2/3 , 



M'(X)\ (12) 



where Ai is the Airy function, Ai'(z) = dAi(z)/dz, and 
X is defined as 



X 



2(V - E) 



(h\vv\) 



2/3 ' 



(13) 



Berry [11] has obtained the first term in Eq. (12) by fo- 
cusing his attention on the coordinate space region close 
to the energy contour V(q) = E. The second term in 
Eq. (12), however, is of the same order of Ti and must be 
included in an asymptotic expansion in k: 

(1) In the energetically allowed region we get (X <C 0): 



A, 



Zirn I 



x cos —\X 



,3/2 
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(14) 



The first term is the Thomas-Fermi spatial proba- 
bility density. 

(2) Outside of the energetically allowed region we get 
(X » 0): 

A^ (q , E ) = ^-^ x -m exp ( _ 2 x 3/2 + ^ 

(15) 
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B. Periodic Orbit contribution 

The oscillatory part of the spectral probability density 
A osc (q, E) is given by 

1 r°° 

A osc (q,£) = — Re / dt(q|e-' Ht / B |q)e iBt / fi . (16) 

A semiclassical expression is obtained by inserting the 
Van Vleck-Gutzwiller propagator in Eq. (16) and evalu- 
ating the ensuing time integral by stationary phase. We 
get [1,2]: 



2 



A osc ( q ,i;) 



(17) 



The sum is over all classical trajectories of energy E that 
start at q and end at q in a nonzero time. The action 
S{<\, q, E) is given by 



S(q",q',£)= f 1 P(q)rfq 



(18) 



with q" = q' 
determinant 



q. The amplitude term involves the 



d 2 S 
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£>(q,£) = (-l)Det| 
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The index /u(q) is equal to the number of times the deter- 
minant D becomes singular along the classical trajectory. 
To make this sum more tractable a smoothing over the 
spatial coordinates q and over the energy E is done: 

a. Coordinate smoothing. We smooth Eq. (17) over 
the spatial coordinates q and evaluate this smoothing 
by stationary phase. The stationary phase condition de- 
mands that the initial and final momentum of the classi- 
cal trajectory be equal. Periodic orbits and orbits in their 
immediate vicinity therefore provide the dominant con- 
tributions to the coordinate smoothing integral. In the 
following we shall assume that this coordinate smoothing 
has been done and was sufficient to wash out the contri- 
bution of orbits not in the vicinity of periodic orbits. For 
a given periodic orbit it is then convenient to use a special 
coordinate system which moves along the trajectory [1]: 
the qi coordinate is chosen along the trajectory and the 
Q2 coordinate is chosen perpendicular to it. The oscillat- 
ing part of the spectral probability density then becomes 
[2]: 

A ose (q,£) = 2 ,, /9 ]P- 1 1 



(2ttS) 3/2 ^ Isil 

(\ lci ltr(M)-2 2 . 
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(20) 



where M is the 2x2 submatrix of the monodromy matrix 
involving coordinates g 2 and p 2 , and mi 2 is one of its off- 
diagonal elements. The index /u(q) is given by [12] 



Mq) 



if 
if 



sign(mi 2 ) 
sign(mi 2 ) 



sign(tr(.M)) 
-sign(tr(.M)) 



where \i m is the Maslov index of the orbit [13] and 
sign(a) = a/ 1 a | . The multiplicity of the periodic orbit 
is denoted by a and is dependent on the type of orbit: 

(1) cr = 2 for a libration (self-retracing orbit), 



(2) (7=1 for a rotation, where the time reversed im- 
age of a rotation is considered as being a distinct 
periodic orbit. 

b. Energy Smoothing. For a chaotic system, the num- 
ber of periodic orbits with period less than r increases 
exponentially [14]: 



dN 



1 



where dN is the number of periodic orbits with peri- 
ods between r and r + dr and topological entropy a is 
a measure of how chaotic the system is. In order to 
counter this exponential proliferation, we introduce an 
energy smoothing in Eq. (4): 



^E-e)=Y J fe(E - E k )\{q\^ k )\ 2 

k 

= J f e (E-E')A(q,E')dE' 



—Re / d*/ £ (*)(q| 



-iAt/h, \ jEt/n 



q)e^ t / ft (21) 



where f c (E) is a smoothing function and f e (t) is its 
Fourier transform. We assume that f c (E) approaches 
6(E) as e —> 0. The energy smoothing parameter e must 
be chosen large enough so that the periodic orbit theory 
converges for the number of periodic orbits available yet 
small enough so that we can discern details about the 
energy spectrum [2,7,15]. Equation (20) then becomes 



A osc (q, E; e) 



(2ttS) 3/2 ^ 
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The sum in Eq. (22) includes periodic orbits which are 
repeated traversals of more fundamental periodic orbits. 
It is therefore convenient to rewrite this equation as a 
sum over these primitive periodic orbits (ppo): 

r, OO 

A osc (q, E; e) = ——^ £ £ /.(r^) 



(2irh) ^ ppg n= i 



^'-M" 2 )| 



x cos 



\ a z m\ 2 J 



F"4 



(23) 



where the superscript n denotes a quantity which corre- 
sponds to the n-th repetition of a primitive periodic orbit 
and [2]: 



= TIT, 



S (n > = nS, fi (n) = (n- l)fi m + 
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m i2 = \ T^" m i2, tr(M (n y) = A? + A£ 

The oscillating part of the semiclassical density of 
states, d 0SC (E;e), is obtained by taking the trace of the 
oscillating part of the smoothed spectral probability den- 
sity: 

d 0SC (E;e) = J rfq A osc (q, E; e). 



With 



and 



we arrive 



+00 



1 

lihi 



dx 




dg„(a) 



ri 1 ) for rotations 
r( 1 )/2 for librations 



at Gutzwiller's trace formula [1] 
1 00 



ppo n = l 



y/\tr(MW) 



x cos[n( 



S 



:A*m)] (24) 



where r is the period of the underlying primitive orbit. 
The multiplicity of the orbits in the above formula is 1 for 
librations and 1 for rotations (we still consider a rotation 
and its time reversed image as distinct periodic orbits). 



III. PERIODIC ORBITS 

The anharmonic Hamiltonian system given in Eq. (1) 
has been studied extensively by Baranger and Davies [5]. 
They found several periodic orbits using an accelerated 
continuation method called the Monodromy method. At 
energy E = 0.001 the main periodic orbits are the os- 
cillations in the x and y directions. As the energy is 
increased these periodic orbits bifurcate. Hence, by fol- 
lowing these periodic orbits and their bifurcations as the 
energy is increased they were able to find many periodic 
orbit families [5]. However there was no guarantee that 
all periodic orbits were found: islands of periodic orbits 
in the E — r plane can appear. Shooting techniques are 
also impractical when the Hamiltonian system under con- 
sideration is moderately chaotic: an orbit which at some 
initial time is very close to a periodic orbit will be vio- 
lently separated from that orbit by the latter's unstable 
manifold. 

The method presented below guarantees that all pe- 
riodic orbits with period less than a chosen maximum 
are found [16]. It complements other methods which 
are highly successful in finding periodic orbits when 
the Hamiltonian system is regular or slightly irregular, 
but break down when the system becomes moderately 
chaotic. We first introduce this method and then apply 
it to the Nelson Hamiltonian. 



A. The Surface of Section Method 

Consider a two-dimensional autonomous Hamiltonian 
system with Hamiltonian H(q\,q2,pi,P2)- Classical tra- 
jectories satisfy Hamilton's equations of motion: 



dH 

dpi ' 



Pi 



dH 
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and although the associated phase space is four- 
dimensional, these trajectories are restricted to the 
three-dimensional energy shell. A suitably chosen two- 
dimensional surface transversal to the phase space flow 
of trajectories can therefore provide a great deal of in- 
formation. Such a surface is called a Poincare sur- 
face of section [17]. The successive piercing of this 
surface by a trajectory produces a sequence of points 
{. . . , X-i, xo, x\, X2, ■ ■ ■} which defines the Poincare map 
P : Xi 1— ► Xi + \. The problem of finding periodic orbits is 
then reduced to one of finding fixed points of the Poincare 
map P and of its repetitions P k . 

Consider a simple closed curve C on the surface of 
section and denote by C the curve resulting from the 
Poincare map P of C. This new curve C will be greatly 
distorted by the stable and unstable manifolds of the pe- 
riodic orbits enclosed by C: C will be stretched along the 
unstable manifolds and correspondingly compressed in 
the direction of the stable manifolds. A necessary condi- 
tion for the existence of a fixed point of P to be enclosed 
by C is that there must be an overlap between the area 
enclosed by C and the area enclosed by C . If there is 
no such overlap there cannot be a fixed point of P in 
the area enclosed by C or by C . Similarly we obtain the 
curve C" as the map _P _1 of C. Then a necessary, but not 
sufficient, condition for the existence of a fixed point in 
C is that there must be a common overlap between the 
areas enclosed by C, C , and C" . 

For a chaotic system this common area will be much 
smaller then the original search area enclosed by C. The 
Poincare-Cartan integral invariant [17] imposes the fol- 
lowing restriction on the projected areas enclosed by the 
curves C, C , and C": 



For a periodic orbit with Lyapunov exponent A and pe- 
riod r the reduction in the original search area is of order 
exp(— 2At). This method therefore works best for a sub- 
stantially chaotic system, which is exactly when conven- 
tional periodic orbit finders fail. A shooting algorithm 
applied to the much reduced search area can then be 
used to zoom in on the fixed point, if it exists. 

Hence by choosing an appropriate set of curves Ci 
whose enclosed areas cover the energetically allowed re- 
gion of the surface of section, we can find all the fixed 
points of P. Extension of this method to find the fixed 
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FIG. 1. The perimeter of a small square, centered at 
y = 0.21 and p y = —0.22, is propagated forward and back- 
ward in time by the Poincare map P and the map P _1 , re- 
spectively. The energetically allowed region is delimited by 
the dashed curve, given by Eq. (25). 

points of P k is immediate. We apply the above ideas to 
the Nelson Hamiltonian, given in Eq. (1). 

B. Periodic Orbits of the Nelson Hamiltonian 

We used the Poincare surface of section parameterized 
by y and p y , with x = and p x > given by energy 
conservation: , 




y 

FIG. 2. An enlargement of the region of interest of Fig 1. 
The small boxes denote the map P of some points on the 
square and the triangles denote the map P -1 of some other 
points on the square. The reduction in search area is by a 
factor of ~ 10 6 



Apart from the vertical libration along the y axis, all 
orbits intersect this surface at distinct points. The energy 
was chosen to be E = 0.8. At these energies, the classical 
phase space is almost totally chaotic. We subdivided the 
energetically allowed region inscribed by the ellipse 

e = § + v 2 (25) 

into small squares of side 0.02. To efficiently propagate a 
side of such a square we subdivided it into unequidistant 
points, with a minimum separation of 10 -7 and a maxi- 
mum separation of 10 -2 . The distance between consecu- 
tive points was adjusted according to the stability of the 
orbits, how close the resulting propagated points were to 
the initial square, and how co-linear was the propagation 
of three consecutive points. 




-4-3-2-10 1 2 3 4 
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FIG. 3. The asymmetric rotation, orbit 43 in Table I, 
whose fixed point was found from Fig. (1) and (2). 

In Fig. 1 we show such a small square and its prop- 
agation forward and backward in time by the Poincare 
map P and the map -P -1 , respectively. In Fig. 2 we en- 
larged the region of interest. Note how the square as 
been so distorted that it effectively becomes a thin line 
and how much smaller the new search area is compared 
to the area of the initial square. In Fig. 3 we show the 
unstable asymmetric rotation (orbit 43 in table 1) whose 
fixed point was found with the help of the two previous 
figures. Table 1 lists the first 80 periodic orbits which 
were obtained by finding the fixed points of the Poincare 
map P and the map P 2 . Besides its period we list what 
type of orbit it is, the trace of its monodromy matrix, 
and a point of intersection of the orbit with the surface 
of section. Although the maps P and P 2 do not have any 
stable fixed points, we believe that there exists an n > 2 
such that the map P n will have stable fixed points. These 
stable fixed points would correspond to orbits with rela- 
tively long periods [18]. 
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TABLE I. The first 80 periodic orbits of the Neison Hamiitonian at energy E = 0.8. We fist whether an orbit is a symmetric 
iibration (SL), asymmetric iibration (AL), symmetric rotation (SR) or asymmetric rotation (AR). We aiso give its period, the 
trace of its Monodromy matrix and a point of intersection the orbit makes with the surface of section. 



no. 


type 


period 


ti(M) 




y 


py 


no. 


type 


period 


ti(M) 




y 


py 


1 


SL 


4.44 


-10.6 


-0 


.40000 


1.13137 


41 


AL 


18 


.61 


309.6 


0. 


.17951 


0.48304 


2 


SL 


6.44 


14.7 





.70062 


0.00000 


42 


AL 


18 


.73 


-8735.0 


-0. 


.44322 


1.06915 


3 


AL 


7.14 


-22.3 





.24646 


1.14857 


43 


AR 


18 


.74 


-1147.4 


0. 


.21124 


-0.21652 


4 


SR 


10.51 


-136.6 


-0 


.88225 


0.00000 


44 


SR 


18 


.74 


-11388.1 


-0. 


.87210 


0.00000 


5 


SL 


11.57 


101.4 





.44113 


0.00000 


45 


SR 


18 


.98 


-1396.0 


-0. 


.82779 


0.00000 


6 


AL 


11.60 


373.0 


-0 


.44162 


1.06801 


46 


AR 


19 


.12 


3540.4 


-0. 


.29562 


1.09236 


7 


AR 


12.79 


-263.2 


-0 


.72413 


0.66552 


47 


AR 


19 


.20 


-2628.6 


-0. 


.77290 


0.53836 


8 


AL 


13.66 


525.9 


-0 


.41534 


1.04915 


48 


AL 


19 


.25 


644.9 


0. 


.15879 


0.18951 


9 


AL 


13.66 


175.2 





.35055 


0.21226 


49 


AR 


19 


.32 


-437.2 


-0. 


.46052 


0.90549 


10 


AL 


13.69 


-359.8 





.65708 


0.25149 


50 


SR 


19 


.39 


-20070.8 


-0. 


.39156 


-1.11943 


11 


SL 


14.32 


1066.8 


-0 


.87045 


0.00000 


51 


SL 


19 


.45 


758.1 


0. 


.14911 


0.00000 


12 


AR 


14.70 


-357.9 


-0 


.62007 


0.80885 


52 


AR 


19 


.48 


-3191.2 


-0. 


.85978 


0.27340 


13 


SR 


14.94 


1584.5 


-0 


.32248 


-1.16141 


53 


SR 


19 


.53 


5684.7 


-0. 


.28495 


1.13991 


14 


SR 


15.12 


-600.0 


-0 


.85477 


0.00000 


54 


AR 


19 


.69 


4327.8 


-0. 


.74676 


0.59850 


15 


AL 


15.41 


275.6 





.50670 


0.52290 


55 


AL 


19 


.74 


2428.7 


-0. 


.71776 


0.56791 


16 


AL 


15.42 


605.7 


-0 


.39170 


1.03190 


56 


AL 


19 


.78 


546.3 


-0. 


.32852 


0.98436 


17 


AL 


15.48 


236.5 





.28249 


0.33265 


57 


AL 


19 


.83 


-6575.2 


-0. 


.30616 


1.07891 


18 


SL 


15.79 


344.2 





.26930 


0.00000 


58 


AR 


19 


.89 


4640.7 


-0. 


.83349 


0.39428 


19 


AL 


16.04 


-4505.3 


-0 


.36982 


1.11927 


59 


AL 


20 


.02 


318.2 


0. 


.13613 


0.53865 


20 


AL 


16.28 


458.2 





.15526 


1.11129 


60 


AL 


20 


.23 


3518.1 


-0. 


.81705 


0.17590 


21 


AL 


16.28 


-1123.5 





.21082 


-1.20070 


61 


AR 


20 


.28 


-1338.1 


-0. 


.73856 


0.47795 


22 


AR 


16.38 


-416.8 


-0 


.55354 


0.86450 


62 


AR 


20 


.65 


-404.6 


-0. 


.42170 


0.91194 


23 


AL 


16.38 


1673.9 


-0 


.81428 


0.36278 


63 


AR 


20 


.68 


-1815.2 


-0. 


.80805 


0.14789 


24 


SR 


16.93 


-1619.0 


-0 


.88025 


0.00000 


64 


AL 


20 


.75 


765.5 


0. 


.11765 


0.24672 


25 


AL 


16.99 


626.8 


-0 


.36992 


1.01575 


65 


AL 


21 


.03 


464.3 


-0. 


.30751 


0.96806 


26 


AR 


17.05 


-900.6 


-0 


.81852 


0.26813 


66 


AL 


21 


.09 


990.0 


0. 


.10225 


0.07501 


27 


AL 


17.11 


281.9 





.22726 


0.41700 


67 


AL 


21 


.19 


2613.6 


-0. 


.68410 


0.60906 


28 


AR 


17.22 


2812.4 


-0 


.31993 


1.11887 


68 


AL 


21 


.35 


306.5 


0. 


.09494 


0.58821 


29 


AL 


17.54 


154.3 





.34423 


0.77392 


69 


AR 


21 


.70 


-1475.9 


-0. 


.70618 


0.52906 


30 


AL 


17.61 


503.3 





.20776 


0.11336 


70 


AL 


21 


.82 


4093.1 


-0. 


.78893 


0.27707 


31 


AR 


17.65 


3552.7 


-0 


.85574 


0.27123 


71 


AR 


21 


.91 


-346.6 


-0. 


.38399 


0.91320 


32 


AR 


17.91 


-442.3 


-0 


.50311 


0.89155 


72 


SL 


22 


.00 


4599.7 


-0. 


.81335 


0.00000 


33 


AL 


17.93 


108.7 





.17521 


1.02050 


73 


AL 


22. 


.15 


863.0 


0. 


.08168 


0.29270 


34 


SR 


17.95 


970.8 





.41123 


0.00000 


74 


AL 


22. 


.22 


364.6 


-0. 


.28530 


0.95055 


35 


AL 


18.08 


-6074.1 


-0 


.33694 


1.09863 


75 


AR 


22. 


.22 


-2152.4 


-0. 


.78385 


0.24099 


36 


AL 


18.14 


5059.3 


-0 


.44865 


1.07300 


76 


SR 


22. 


.39 


-2381.0 


-0. 


.80242 


0.00000 


37 


AL 


18.15 


2123.6 


-0 


.75958 


0.49887 


77 


AL 


22. 


.55 


2697.8 


-0. 


.65554 


0.63591 


38 


AL 


18.37 


-1752.4 





.20504 


-1.20899 


78 


AL 


22. 


.60 


1192.8 


0. 


.06291 


0.13077 


39 


AL 


18.44 


603.2 


-0 


.34910 


1.00008 


79 


AL 


22. 


.62 


272.6 


0. 


.05395 


0.63485 


40 


SL 


18.45 


2716.4 


-0 


.84063 


0.00000 


80 


SL 


22. 


.74 


1304.0 


0. 


.05667 


0.00000 



SEMICLASSICAL RESULTS 

The calculation of the exact quantal wave functions 
and energies were obtained by matrix diagonalization. 
We used the basis <f>m\%)<f>nf\y — % 2 /2), where (f)^ and 
are harmonic oscillator wave functions appropriate 
for low energies. The matrix elements of the Nelson 
Hamiitonian in this basis are given in Appendix A. We 
set k = 0.05 and used 240 oscillator states for the x di- 
rection and 26 for the y direction. The maximal number 
of basis states in the x direction and the y direction were 
chosen in conjunction with the value of k in order to 



get the best possible convergence of the eigenvalues up 
E ~ 0.8. At these energies the classical phase space is 
for all practical purposes totally chaotic. The accu- 
racy of this basis was further tested against the average 
semiclassical spectral staircase N avg (E), which is given 
by (see Appendix B) 

(26) 

The first term is the usual Thomas-Fermi expression for 
the staircase. The inclusion of corrections to the Thomas- 
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275 
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FIG. 4. Comparison of the average semiclassical staircase, 
and the one obtained by diagonalization. The dashed line is 
the Thomas- Fermi staircase. 

Fermi staircase is essential to achieve a proper fit with 
the exact spectral staircase. This is shown in Fig. 4 for 
the energy interval 0.7 < E < 0.8. The agreement be- 
tween the semiclassical prediction for the average stair- 
case N avg (E) and the exact staircase is excellent. In 
Fig. 5 we show the oscillating part of the staircase func- 
tion N 0SC (E) = N(E) - N avg (E). The sudden dip in 
Nosc(E) around E ~ 0.83 indicates the failure of the 
diagonalization for higher energies. 




FIG. 5. The oscillatory part of the spectral staircase, 
Nosc(E) = N(E) - Navg(E). The sudden dip indicates the 
failure of the diagonalization for E ~ 0.83. 




0.6 0.65 0.7 0.75 0.8 



C. Density of States 

Gutzwiller's trace formula, Eq. (24), provides the 
means of calculating semiclassically the oscillating part 
of the density of states d 0SC (E; e). We chose the following 
form for the energy smoothing: 



ME) 



1 



-E 2 /2e 



FIG. 6. The oscillating part of the smoothed energy level 
density, d osc (E; e), for three different energy smoothing pa- 
rameters e. The dashed lines represent the semiclassical re- 
sults while the solid lines display the exact results. The aver- 
age energy level density at E = 0.65 is ~ 580. Top: e = 0.02, 
and the first 2 periodic orbits in Table I contribute noticeably 
to the sum in Gutzwiller's trace formula. Middle: e = 0.01, 
and the first 5 periodic orbits contribute. Bottom: e = 0.005, 
and the first 60 periodic orbits contribute. 

The damping of the long orbits in Gutzwiller's trace for- 
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mula is then given by 

/ e (r) = e-^ 2 / 2s2 . 

We applied Gutz wilier 's trace formula with different 
values of the energy smoothing parameter e. In Fig. 6 we 
show our semiclassical results (dashed lines) and compare 
them with the exact calculations (solid lines). We took 
e = 0.02 in Fig. 6a. Only the first 2 periodic orbits listed 
in Table 1 contribute noticeably to d 0SC (E;e). Since the 
average semiclassical density of states at E = 0.65 is 
~ 580, the maximum oscillations in the density of states 
is about 0.4% of the average value. The agreement be- 
tween the exact and semiclassical results is excellent up 
to energy E = 0.75. It is the failure of the numeri- 
cal diagonalization for energies above E = 0.83 that is 
causing the exact results for d 0SC (E;e) to diverge after 
E = 0.75. Note that the number of states in an en- 
ergy range [E — e, E + e] is about 24 for e = 0.02. Thus 
Gutzwiller's trace formula only describes the gross fea- 
tures of the clustering of the eigenvalues for this value of 
e, but does so extremely well. 




-2-10 1 2 



X 

FIG. 7. The oscillating part , A osc , of the energy and spa- 
ce smoothed coordinate space density. Energy smoothing pa- 
rameter e = 0.01. Top: exact quant al. Bottom: Eq. (23) 
with 5 orbits. Solid contours are positive, dashed are nega- 
tive. The contour spacing is 0.5. The Thomas-Fermi density 
is 63.7. The energy is E = 0.692. 

In Fig. 6b we took e = 0.01. In this case the first 5 
periodic orbits contribute significantly and the maximum 
oscillations are about 3% of the average value and there 



are about 12 states within an energy range [E — e, E + e]. 
The agreement is still very good. However, there are 
some small discrepancies and the inclusion of more peri- 
odic orbits do not help. This is also the case in Fig. 6c 
where e = 0.005. Only the first 60 periodic orbits con- 
tribute and now the maximum oscillations are 13% of the 
average value and about 6 states are in an energy range 
[E — e, E + e\. Since Gutzwiller's trace formula is the 
first term in an expansion in powers of Ti of d 0SC (E; e), we 
believe that higher order corrections are needed to get 
better agreement between the exact and semiclassical re- 
sults when e becomes smaller. 



D. Spectral Probability Density 

Bogomolny's formula, Eq. (23), provides us with the 
means of calculating semiclassically the oscillatory part of 
the spectral probability density A osc (q, E; e). In Figs. 7, 
8, and 9 we show our semiclassical results (bottom half 
of figures) and compare them with the exact calculations 
(top half of figures). The solid contours indicate when 
A osc (q, E; e) is positive and the dashed contours when 
A osc (q, E; e) is negative. The contour spacings is 0.5 in 
all Figs. 7, 8, and 9. 

The energy smoothing parameter was set to e = 0.01 
and we found that only the first 5 periodic orbits con- 
tribute significantly. We also performed a gaussian 
smoothing over the coordinate space with 

J_ e -(««+«J)/» a 

We took b = 0.2, which is sufficient to wash out the con- 
tribution of most non-periodic orbits in the central re- 
gion of the potential. However this value of the smooth- 
ing parameter b is too small for the "arms" region of the 
potential. The nodal structure seen up the arms of the 
potential can easily be described using WKB expressions 
for the wavefunctions since the classical motion up the 
arms is adiabatically separable. We therefore concen- 
trated our attention to the central part of the potential. 
The chosen value of & is a compromise between choosing b 
small enough to see details in the central region and large 
enough to eliminate the contribution of non-periodic or- 
bits to A osc (q, E; e). It is interesting to note that the 
gaussian coordinate smoothing of the spectral probabil- 
ity density is equivalent to the coordinate space projec- 
tion of the coherent state representation of the spectral 
operator [19]. 

The average part of the coordinate smoothed and en- 
ergy smoothed spectral probability density was calcu- 
lated using both Eq. (12) and Eq. (CI) from Appendix 
C. It was found that Eq. (CI) sufficed for the present 
accuracy of the calculations and was much faster. The 
Thomas-Fermi density is 63.7. The maximum oscilla- 
tions in A osc (q, E; e) were around 10 or about 15% of 
A avg (q,E;e). 
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-2-10 1 2 

X 

FIG. 8. Same as Fig. 7 for E = 0.700. 




X 

FIG. 9. Same as Fig. 7 for E = 0.708. 




0.5 0.6 0.7 0.8 

E 

FIG. 10. The exact scar strength A (full line) and the 
semiclassical one calculated with the first 5 orbits in Table 
I (dashed), as a function of energy, at three points of coor- 
dinate space. Bottom: x = 0, y = —0.65, emphasizing orbit 
no. 1. Middle: x = 0.8, y = 0, emphasizing orbit no. 2. Top: 
x = 0.75, y = 0.75, emphasizing orbit no. 3. 

The qualitative agreement between the exact quantal 
and the semiclassical pictures is very good. The main 
trends are well represented in all of our comparisons. 
The average density of states at E = 0.70 is ~ 626. 
Hence there is a shift by about 5 eigenstates between 
the calculations presented in Fig. 7 and Fig. 8 and simi- 
larly between Fig. 8 and Fig. 9. Since the chosen energy 
smoothing parameter (e = 0.01) is relatively large, we 
were surprised that the drastic changes in A osc (q, E;e) 
are as well represented by the semiclassical calculations. 
In Fig. 10 we show A osc (q, E; e) for three points in co- 
ordinate space for 0.5 < E < 0.8. The three points are 
chosen so as to emphasize separately the contribution of 
the three lowest periodic orbits. Again the agreement is 
very reasonable considering the approximations. 

IV. CONCLUSION 

We showed that periodic orbits play a central role in 
the semiclassical description of a physical system. The 
method we used to find the periodic orbits guaranteed 
that we had all the periodic orbits up to a chosen maxi- 
mum period. This allowed us to apply Gutzwiller's trace 
formula and Bogomolny's formula to a generic Hamilto- 
nian system. 
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In our calculations of the oscillating part of the den- 
sity of states, we found that when the energy smooth- 
ing parameter was large enough Gutzwiller's trace for- 
mula worked extremely well. Not only did we get a very 
accurate semiclassical description of the oscillations in 
the density of states, we did so with little numerical ef- 
fort since only a few periodic orbits were needed. Also, 
exact quantum calculations become numerically infeasi- 
ble for higher energies. Gutzwiller's trace formula then 
is invaluable in describing the coarse oscillations in the 
density of states. However when we reduced the energy 
smoothing parameter so as to obtain more detailed in- 
formation about the spectrum, we found that there was 
a discrepancy between the amplitudes of the semiclassi- 
cal and exact quantal calculations. These discrepancies 
could not be reduced by including more periodic orbits 
in Gutzwiller's trace formula. We noted that because the 
amplitude of the oscillations relative to the average den- 
sity of states was not small, asymptotic corrections in k 
to Gutzwiller's trace formula are needed. 

Similar calculations of the oscillatory part of the spec- 
tral probability density using Bogomolny's formula were 
done. We found that there was very good qualitative 
agreement between the exact quantal and the semiclas- 
sical calculations over a wide energy range. Since the 
maximum oscillations in A osc (q, E; e) were around 10 or 
about 15% of A a „ 3 (q, E; e), better quantitative agree- 
ment will only be possible if asymptotic corrections in Ti 
are included. 
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and then to make the change of variables from {x, y} to 
{6,6}: 

~ + 6 2 )£* - (26^ + l^Jnonumber (Al) 



^ t 2 _i_ f 2 



^(6,6) = ^(6,6) 



The matrix elements of the Hamiltonian in the harmonic 
oscillator basis with the variables £1 and ^2 are obtained 
by making the identification 

h d h d 

1 dix 1 d£ 2 



and rewriting Eq. (A2) as 



f + (l+€l 2 )f -&Pl-f)P2+ftf+6 



= Ei> 

With wi = and w 2 = \/2, 



Tiwi 



Pi = - ai) ^ 1 = ^2^^ 1 + a ^ 



V2 = y -rr l { a 2 - a 2) 6 



— («2 + «2) 



the desired matrix elements of H in the basis \ri1ri2) 
l n 2)f 2 are obtained: 

(ni n 2 \H\n 1 n 2 ) = TiL0i{rii + ^) + ^w 2 (n 2 + ^) 



(A2) 



APPENDIX A: DIAGONALIZATION 

The eigenvalues and the corresponding eigenstates of 
the Nelson Hamiltonian H were obtained by diagonal- 
ization in a displaced harmonic oscillator basis with the 
spatial coordinates £1 and £ 2 : 



6 



6 



The coordinate transformation from {x,y} to {£i,£ 2 } 
is non-orthogonal. A simple way to make this coordi- 
nate transformation is to write Schrodinger's equation in 
Cartesian coordinates, 



' 2 



dx' 2 



dy' 2 



} + F + («-t) 



ip(x,y) 



Eip(x,y) 



(nin 2 \H\ni 7i 2 +2) 
(rii n 2 \H\rii+2 n 2 ) 



■^-(ni + i)V(n 2 + l)(n 2 + 2) 



Tl 2 L0 2/ 1, 



4 Wl " (n2+ 2 )A/(ni + 1)( " 1+2) 



117121^1711+2 712+1)= fy^V("2 + 1)("1 + l)(7ll + 2) 



U n 2 \H\ ni +2 712-1) = --y — \/n 2 { ni + l)(7li + 2) 



(tIi71 2 |#K+2 7l 2 +2) 
(Til 712 1^1^1+2 7l 2 -2) 



Tl 2 LU 2 

' 8wi 

Tl 2 LU 2 

' 8wi 



V(7l2 + l)(7l2+2)(7li+l)(7li+2) 

\/n 2 {n 2 -l)(ni + l)(ni + 2) 
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APPENDIX B: AVERAGE DENSITY OF STATES 



which simplifies to 



The Thomas-Fermi approximation for the staircase 
N(E) = '^2 i 6(E — E{) is the number of Planck cells con- 
tained in the phase space volume enclosed by the energy 
shell. For the Nelson Hamiltonian, 



N aV g(E) 



1 



{2-why 



A* 2 



rfqrfp 6 IE 



Px 



P£_ 
2 



r~ 2 / \2 

77* -\V--z-) 



The set of transformations 



12 x 2 

x = \i-zi, y--^ 



*2, P: 



X = V2 



Z 3 , Py = \/2z 4 



(Bl) 



Taking the inverse Fourier transform of Z(p)//3 we get 
IE 2 ' * 2 



Navg{E) 



h 2 



Ti 2 + 2 2 i , 4l 
1 -l2[V + ^] + °^ 



Finally the derivative of N avg (E) with respect to E 
gives the average level density d avg (E): 



APPENDIX C: AVERAGE SPECTRAL 
PROBABILITY DENSITY 



N avg (E) 



A* (2TTh) 2 

J d 4 z e(E- Zl 2 - z 2 2 - z 3 2 - z 4 2 ) 



The above integral is the volume of a four dimensional 
sphere of radius \f~E~ . Therefore 

1 E 2 

Navg{E) 



2^l n 2 

To obtain corrections to N avg (E), we define Z(/3) as the 
trace of the evolution operator. From Eqs. (8) and (11) 
we get for the Nelson Hamiltonian, up to 0{Ti 3 ): 



z(p) 



(2irhY 



dqdpe 



x(i- y A 



with 



A\ SS (x, y,f)) = \p 2 {(^ + 2) + 2x 2 -2(y-^-)} 

-{x 2 ^ 2 -4n(y- y) +*n(y- 

The set of transformations (Bl) gives 



3 



x 
~2 



z(p) 



1 



a* (2Trhy 



dz\ dz2 dzs dz4 e 



+Z2 +Z3 +Z4 



P 3 r2zi s 

3 v ii 



^ 2 + W)+W] ] } 



We found during our numerical work on the spatial 
wavefunctions that coordinate and energy smoothing of 
the Thomas-Fermi spatial density was sufficient to re- 
produce the average spectral probability density. The 
Thomas-Fermi spatial density is 



A TF (q,E) 



2^ lf V(<l)<E 



{ otherwise 
With a gaussian smoothing in energy, we get 
1 A 00 

Aa„,(q,£)=— 2 / UE'-E)dE' 

/•oo 

■. / exp \-{E' - Ef/2e 2 \dE' 

Jv(a) 1 J 

(CI) 



1 1 
2ttS 2 V2^? 



where erfc is the complementary error function. The co- 
ordinate smoothing was then done numerically. 
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